Black brane solutions of Einstein-Maxwell-scalar 
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Abstract 

We investigate the global properties of black brane solutions of a three-parameter Einstein- 
Maxwell model nonminimally coupled to a scalar with exponential potential. The black 
brane solutions of this model have recently been investigated because of their relevance 
for holography and for the AdS/condensed matter correspondence. We classify all the 
possible regular solutions and show that their asymptotic behavior either breaks hyper- 
scaling invariance or is that of a domain wall. We also write down some exact solutions in 
! Schwarzschild coordinates. 
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1. The model 

The Einstein-Maxwell model nonminimally coupled to a scalar field with Liouville 
potential (EML) has been studied in several contexts [1-5]. Models of this kind arise for 
example as low-energy limits of supergravity and superstring theories [6]. They are espe- 
cially relevant in the context of holography [3-4] , since they are dual to models exhibiting 
a breaking of the conformal symmetry, leading to hyperscaling violations [7] . 

Because of the nonminimal coupling, the standard no-hair theorems [8] do not hold 
for these models, and their black hole solutions have unusual asymptotics and support 
nontrivial scalar fields. In particular, in the case of spherical symmetry, the possible 
black hole solutions have been classified in ref. [1] , where it was shown that all spherically 
symmetric solutions have non-standard asymptotic behavior. Because the potential for 
the scalar field does not admit local extrema, these models do not allow for AdS vacua. 
Hence, also the Reissner-Nordstrom-AdS black hole solution is forbidden. It should be 
noticed that instead AdS2 x S 2 solutions could be present. 

Particularly interesting are also the planar solutions, in particular black branes. From 
an analytical point of view, planar solutions have the advantage of being easier to find in 
closed form than spherically symmetric ones, since they require one less integration, but 
are especially interesting because of their relevance for holography and the AdS/condensed 
matter correspondence [3,4]. 

Although EML models do not allow for AdS vacua, they still could admit asymptotic 
solutions preserving scaling isometries. These solutions would be asymptotic to Lifshitz 
spacetimes [9,10], i.e. spacetimes presenting a scale isometry under which timelike and 
spacelike coordinates transform with different exponents. This property may be relevant 
for the holographic description of quantum phase transitions. 

In the general case, the scale isometry is broken and the metrics only transform 
covariantly under scale transformations [5]. To this class belong in particular solutions 
for which Poincare invariance of the brane is preserved. In the literature they are often 
referred to as domain wall solutions [11,4,12-14], in analogy with the domain wall solutions 
of supergravity theories. In a EML context they were studied in [4]. In a wide range of 
the parameters, domain wall solutions are conformal to AdS and allow for an holographic 
interpretation [11,4]. 

It has been shown that the holographic interpretation of Einstein-Maxwell-scalar the- 
ories gives rise to a very rich phenomenology in the dual QFT. This includes phase transi- 
tions triggered by scalar condensates and non-trivial transport properties of the dual field 
theory [2-5,10,12,15]. From the holographic point of view the difficulty connected with 
the absence of an AdS vacuum, and hence of an ultraviolet fixed point in the dual QFT 
can easily be circumvented. In fact EML models can be considered as the near-horizon, 
near extremal description of Einstein-Maxwell-dilaton gravity whose scalar field potential 
allows for a local maximum, hence for an ultraviolet fixed point [16]. 

Very recently, it has been realized that the solutions of EML theory represent the 
main example of scale covariant models that lead to hyperscaling violation in the dual field 
theory [7] . They are therefore a very promising framework for the holographic description 
of hyperscaling violation in condensed matter critical systems (e.g. Ising models [17]). 
Moreover they have been used for the description of Fermi surfaces and for the related 



2 



area-law violation of entanglement entropy [7,18]. 

Although planar solutions of the EML model have been already investigated and 
several partial results have been obtained [1-4], in this paper we systematize the previous 
findings and obtain general results for all values of the parameters of the model. In 
particular, we study the phase space of the dynamical system associated with the planar 
solutions of the model in four dimensions, using the same approach adopted in ref. [1] for 
the spherically symmetric case, and classify all the possible black brane solutions admitting 
regular horizon, according to their asymptotic behavior. It turns out that all solutions have 
nonstandard asymptotics, of domain wall or hyperscaling-violating form, but not of anti-de 
Sitter or Lifshitz form. Moreover, for given values of the parameters of the action, solutions 
can exist that behave differently asymptotically. 

We also obtain some exact solutions for special values of the parameters. Some of 
these solutions have already been obtained in [2], using a different parametrization of the 
metric, that gives rise to complicated field equations, and less clear interpretation. 

Since the main purpose of this paper is the classification of the black brane solutions 
of the EML theory with regular horizon and asymptotic region, we shall not investigate 
other possible solutions (namely solutions containing naked singularities or cosmological 
horizons). We also do not discuss the physical applications of the model, for example to 
holography, except for some considerations on thermodynamics. A thorough investigation 
of these topics can be found in [3-5] . 

2. Action and field equations 

We consider the action 



where F^ v is the Maxwell field, is a scalar field and g and h are two real parameters. 
Some special cases are well known: for h — > oo one obtains the GHS model [19,20], for 
h = the GHS model with a cosmological constant. For (7—7-00 one gets a Liouville model, 
for g = the minimally coupled Einstein-Maxwell theory with a Liouville potential. For 
h = —g = ±1, the action can be derived from string theory [11], for h = —l/g = ±\/3 
from the Kaluza-Klein reduction of a five-dimensional model. 
The field equations read 



We look for electrically charged solutions with planar symmetry. Magnetic solutions 
can be obtained by duality. More precisely, electrically charged solutions of the model with 
parameter g are magnetically charged solutions of the model with parameter —g [2]. 




(2.1) 




-2hcf> 



(2.2) 
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It is useful to parametrize the metric as* 

ds 2 = -e 2v dt 2 + e 2u+Ap de + e 2p (dx 2 + dy 2 ), (2.3) 

where v, p and 4> are functions of the radial coordinate £. 
In these coordinates, the Maxwell equation is solved by 

F H = e 2 ^+^Q, (2.4) 

with Q the electric charge, and the remaining field equations can then be put in the form 

v" = Ae 277 + Q 2 e 2 *, (2.5) 

p" = \e 2 ^ - Q 2 e 2 *, (2.6) 

0" = h\e 2ri + gQ 2 e 2x , (2.7) 

p' 2 + 2pV - <P' 2 - Ae 2 " + Q 2 e 2 * = 0, (2.8) 

where 

T] = v + 2p - h(j), x = v + g ( j ) . (2.9) 

Writing the equations (2.5), (2.7) and (2.8) in terms of the independent variables 1], x an d 
p, one obtains 



v " = (3 - h 2 )Xe 271 - (1 + gh)Q 2 e 2x , (2.10) 
X " = (1 + gh)Xe 2 ^ + (1 + s 2 )Q 2 e 2 *, (2.11) 
ap' 2 + X' 2 + V' 2 ~ 2r/x' - + Wx* + (9 + h) 2 (Xe 2 ^ - Q 2 e 2 *) = 0, (2.12) 



where 

a = 4 + 3g 2 + 2gh-h 2 , 71 = 2 + g 2 + gh, 72 = 2 - h 2 - gh. (2.13) 

The case a = is singular, since (2.10) and (2.11) are no longer independent. 
The system (2.5)-(2.8) is invariant under the shift 

v^tv-ge, 4> ->• + e, (2.14) 

with constant parameter e. This invariance, which is not present in the spherically sym- 
metric case, facilitates the solution of the system. In fact, it implies the existence of a 
conserved quantity, that can be obtained combining (2.6), (2.10) and (2.11). Indeed, 



p = , (2.15) 



a 



* This parametrization was originally introduced in [20] and then used in several papers 
for discussing the phase space of solutions of more complicated models [21,1]. 



and hence 

lW ~ I2X' ~ ap' = b, (2.16) 
with b an integration constant. Substituting in (2.12) one gets 

b 2 



(g + h) 2 



+ (1 + g 2 ) V ' 2 + 2(1 + gh)rj' X ' - (3 - h 2 ) X ' 2 - a(Ae 2 " - Q 2 e 2 *) = 0. (2.17) 



Eqs. (2.10) and (2.11) form a dynamical system for the variables r\ and x, subject to the 
constraint (2.17). 

3. Exact solutions 

In a few special cases the field equations can be solved exactly. These solutions are 
useful for the discussion of the general case. 

A. Neutral solutions, Q = 

A simple case in which the field equations can easily be solved exactly is when the 
electric charge vanishes, Q = 0. This is a minimally coupled gravity-scalar model with 
exponential potential. The standard no hair theorems [8] do not apply because the solutions 
are not spherically symmetric, and hence nontrivial solutions are possible. 

The field equations (2.5)-(2.8) reduce to 

v" = p" = ^" = \e 2 \ (3.1) 
p' 2 + Iv'p' - c/>' 2 - Ae 2 " = 0. (3.2) 

In addition to (2.14), the system now enjoys a further symmetry, for v — > v+k, p — > p+h 
(f) — > (f) + 2k, with constant parameter that permits to completely integrate the system. 
From (3.1) it follows that 

, V + b , i]' + d , hrf + c 

V = 3^' P = 3^7?' = 3^' (3 - 3) 

with integration constants 6, c, d, and c = (b + 2d)/h. Solving (2.10) with Q = yields 

,' 2 = (3-ft 2 )Ae^ + « 2 , >^= (3 _^^ r M 

Substituting (3.3) and (3.4) into (3.2), one gets a constraint between the parameters of 
the solution, 

h 2 (3 - h 2 )a 2 + 2(h 2 - 2)bd + (h 2 - A)d 2 = b 2 . (3.5) 

The radial function e p ~ e (v+d£)/(3-h ) g oes ^ Q infinity as £ — > 0, and we identify this 
limit with spatial infinity. If the solution possesses a regular horizon, in the opposite limit 
£ — > — oo, e p must be a nonvanishing constant. This request enforces the choice d = —a, 
and hence b = (2 — h 2 )a, c = —ha. 



Integrating (3.3) and substituting the solution (3.4), one can obtain the explicit form 
of the metric functions u, p and <f>. Rather than pursuing with the present coordinates it 
is however convenient, in order to get a more transparent interpretation of the solutions, 
to write the metric in a Schwarzschild gauge, namely, 

ds 2 = -Udt 2 + U^dr 2 + R 2 (dx 2 + dy 2 ), (3.6) 

where U(r) = e 2l/ , R(r) = e p , and F tr = -jjz- In this gauge, the field equations take the 
form 

d 2 R fd£ 2 



dr 2 \ dr 

* ( UR^) = -%e 2 ^ + XR 2 e- 2h t (3.7) 
dr \ dr J R l 

A. ( UR 2 ^-) = g^e 2 ^ + hXR 2 e~ 2h t 
dr \ dr J R 2 

The new coordinate r is related to £ by 

r = J 'e^+ridt « (1 - e ^)-(i+h')/(s-h')^ (3 8) 

where « means modulo a constant factor. It follows that 

e 2 < « 1 - ^-tf-^/a+fc 2 ), e 2r > « (l - ^--(^V^ 2 )) ^(3-/^/(1+/^ (3 . 9) 

with \i a constant. From the solutions obtained above, one can finally write down the 
metric functions in terms of r and the new integration constants A, £?, C, /jl, as 

R = e P = Ar 1 ^ 1 ^, U = e 2 " = B ( 1 - ^-O-O/a+ft 2 )) r 2/(i+h 2 ) 

V / (3.10) 

e 2< ^ = C 2 r 2/l/(1+/l2) . 

The constants can be fixed substituting (3.10) in (3.7). It turns out that \i is a free 
parameter, while two of the constants, say A and C, can be set to 1 by rescaling the 
coordinates. Then 

(l + h 2 ) 2 X 



B 



3-h 2 



These solutions have been found in [13] using a different method. 

For positive A and h 2 < 3, the solutions (3.10) describe a one-parameter family of 
black brane solutions with horizon at rn = ^ 1+h ^^ s ~ h ) and domain- wall asymptotics, 
that reduces to anti-de Sitter if h = 0. Solutions exist also for negative A, if h 2 > 3. 
In this case the fi term dominates and an asymptotic region exists only if /i < 0, with 
JJ ~ r -(i-^ )/(i+h )_ However, a naked singularity occurs at r = 0. 
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B. Planar GHS solutions, A = 

Another case in which it is possible to obtain exact solutions is when the potential 
vanishes, A = 0. This is the planar extensions of the well-known GHS solutions [19-20]. 

Also in this case a new symmetry is present that enforces complete integrability. The 
field equations (2.5)- (2.8) reduce to 

v" = -p" = = Q 2 e 2 *, (3.11) 
9 

p' 2 + 2v'p'- ( j ) ' 2 + Q 2 e 2 x=0, (3.12) 

and the system is invariant under (2.14) and u —> v + k, p ^ p, — k. 

From (3.11) follows that 

/ x' +b , -x' + d , gx' + c 

1 + g z 1 + g z 1 + g z 

with c = —hjg, and eq. (2.11) is solved by 

x « = (l + s W + a 2 , Q^= (1 + g 4 2 g 2 °^ )2 . (3.14) 
Moreover, substituting in (3.12) one gets the condition 

b 2 

-(1 + g 2 )a 2 +2bd + d 2 = (3.15) 

As in the previous case, it can be shown that a horizon is present if d = a. It follows 
that b = g 2 a, c = —ga. Defining the coordinate r as in (3.8), with a suitable choice of the 
integration constant one obtains 

r w e 2aC - 1. (3.16) 
Passing to Schwarzschild coordinates (3.6), 

R = e p = Ar 1 ^ 1 ^, U = e 2 ^ = Br~ 2 l^^ (1 + F ), 

(3.17) 

Also in this case, substituting in (3.7), one can easily check that \i is free, A and C can be 
set to 1, and 

B- Q2 
~ 1+1?' 

The properties of the solutions depend on the sign of the free parameter fi, which dic- 
tates the asymptotic behavior. An asymptotic region is present only if p, > 0, with 
U ~ r -0--9 ). i n this case, a horizon is present at r = — p, but a naked singu- 

larity occurs at the origin, so these are not regular black brane solutions. For p < 0, the 
horizon is of cosmo logical type, but a naked singularity is still present. 
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C. Special case x' = v' 

Finally, an exact special solution can be found for generic values of the parameters of 
the action, when x' = V' \ i- e - X = V + log-fT, with constant K. In this case, 

, (g + hfrj + b , (g + h)(4 + g 2 -h 2 ) V ' - 2gb (g + h) 2 V ' + b 

P — > v — — — — , <p —Z — — — , [6. Lb) 

a a(g + h) a{g + h) 



with b an integration constant. Moreover, comparing (2.10) and (2.11), one gets 
and (2.12) becomes 



+ v ' 2 = -\e 2 \ (3.20) 



(g + h)± ' 7i 

with a, 71 and 72 given by (2.13). 

Solving (2.10), with Q 2 e 2x = ^\e 2r >, one then obtains 



71 Q;(l-e 2 <) 2 



(3.21) 



with a an integration constant. Substituting this result in (3.20), one gets b 2 = (g + h) 4 a 2 . 

As discussed before, a horizon is present if e p goes to a constant for £ — >■ —00. This request 
enforces the choice b = —(g + h) 2 a. 

Defining then the radial coordinate r as in (3.8), one obtains 

r ~ (I _ e 2<)-/3/a 5 (3.22) 

where 

P = 4 + (g + h) 2 . (3.23) 
Again, one can write the metric in the form (3.6). One has 

V P J (3.24) 

Q 24> _ C <2 r 4(g+/ l )//3_ 



The constants can be determined substituting in (3.7). In particular, one of them, say 
C, can be set to 1, while it is easy to check that \x is a free parameter. The remaining 
constants result in 

* = B — X. 

y 72 A ct7i 

These solutions have been obtained in different coordinates in [2]. 
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The properties of the solutions depend on the values of the parameters g, h and A. 
For positive A and 0,71,72 > or 0,71,72 < 0, the solutions (3.23) describe a charged 
black brane, with horizon at = and hyperscaling-violating asymptotic behavior, 

that in the special case g = 2/h — h is enhanced to that of a domain wall. If A is negative 
and 71 > 0, 72 < 0, one also has black brane solutions when a > 0. In all other cases no 
asymptotic region is present or naked singularities occur. 

4. Thermodynamics 

Before passing to the study of the global structure of the space of solutions of the 
system (2.10)-(2.12), we calculate the thermodynamical parameters of the exact black 
brane solutions found in the previous section. We omit the A = case, since it contains 
naked singularities. 

A. Q = 0. 

The mass density can be defined as [10] 



m = lim — — Vu 

r— 5>oo SlV 



, j. dR 2 yfjj *^"^ 2 
dr dr 



(4.1) 



where we have subtracted a term corresponding to the solution evaluated on the back- 
ground with rh = 0. For the solutions (3.10) this gives 

m = —r u, (4.2) 

The temperature T can be obtained calculating the periodicity of the Euclidean section, 
as 

1 dU(r h ) lA , 

T=- 4.3) 

An dr 



In our case 



T=± (l + h 2 )ti {1 - h2)/{3 - h2) . (4.4) 

47T 



Finally, the entropy density is simply obtained as 

1 



Hence 



It follows that dm = Tds and m = -^f^ 



s = - A R\r h ). (4.5) 



s = l^/(3-/, 2 ). (4 . 6) 

Ts 
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B. X ' = rf 

Using the previous formulae, one obtains for the solutions (3.23) 



(9 + h) 2 (3 XQ 2 
m= — W /I, 4.7) 

8-7T ay 7172 



and 

T_ A A ..(4+g 2 -2gh-3h 2 )/ a 
A r 1 > 

47T 71 



4 V 72 A 



It is easy to see that dm = Tds and m = 2 ^ 9+h ^ Ts 



5. The dynamical system. 

In the general case, following the methods of ref. [21,1], eqs. (2. 10), (2. 11) and (2.17) 
can be put in the form of a dynamical system, by defining X = x'> Y = rj ', P = Qe x , 
Z = v^AJe 7 ', with 

X' = (l + gh)eZ 2 + (l+g 2 )P 2 , 

Y> = (3 - h 2 )eZ 2 - (1 + gh)P 2 , (5.1) 
Z' = YZ. 

The independent variables are X, Y and Z, while P 2 is defined as 

P 2 eZ 2 (1 + g 2 )Y 2 - (3 - h 2 )X 2 + 2(1 + gh)XY b 2 

a (g + h) 2 a ' ' 

where a is given by (2.13) and e = +1 if A > or e = — 1 if A < 0. For each value of 
the free parameters b and Q, we can now discuss the structure of the phase space of the 
dynamical system. 

The system is invariant for Z — > —Z, so it is sufficient to discuss the Z > portion of 
phase space. Moreover, it is invariant also for X — > —X, Y — > —Y, together with £ — > — £, 
and hence to each critical point with X > corresponds a critical points with X < 0, 
with the direction of the trajectories reversed. Therefore, in the following we list only the 
critical points with positive X. 

The critical points of the dynamical system at finite distance are attained in the limit 
£ — > ±oo, and correspond to finite Schwarzschild radial coordinate r [21]. They lie at 
Z = P = 0, and hence their coordinates X , Y satisfy the equation 

(1 + g 2 )Y 2 - (3 - h 2 )X 2 + 2(1 + gh)X Y = . (5.3) 
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If a > 0, this is the equation of a hyperbola in the Z = plane, while for a < it 
represents an ellipse. In the degenerate case a = 0, the set of critical points is given by 
a pair of straight lines. Another limit case is b = 0, that corresponds to extremal black 
branes: in this case the critical points again lie on a pair of straight lines. In all cases, the 
trajectories starting from the curve of critical points correspond to the limit £ — > — oo and 
those ending on the curve to £ — > oo. 

The characteristic equation for small perturbations around the critical points has 
eigenvalues 0, 2Xq and Yq. Hence, for a given value of 6, each point in the Z = plane 
satisfying (5.3) with Xq > 0, Yq > 0, repels a 2-dimensional bunch of solutions in the full 
3-dimensional phase space. The points with Xq > 0, Yq < act as saddle points. The 
presence of a vanishing eigenvalue is due of course to the fact that there is a continuous 
set of critical points lying on a curve. 

Integrating (2.15) when a^0, one gets the expression for the radial function, 

e p = const, x e (-nv-™-i*)/<*. (5.4) 



In general, e p vanishes as £ — > — oo, except in the special case 71 Yq — 72^0 — b = 0. This 
equation, combined with (5.3), has the only real solution Xq = Yq = b/(g + h) 2 . In this 
case e p — > const as £ — > —00, indicating the presence of a horizon. Therefore, all the 
trajectories starting from the point X = Y = b/(g + h) 2 describe black brane solutions, 
while those starting from the other critical points describe naked singularities. 

To complete the analysis of the phase space one must also investigate the nature of the 
critical points on the surface at infinity, that correspond to the asymptotic region of the 
physical solutions and are attained for a finite value of £. The analysis can be performed 
defining new coordinates u, y, and z such that the surface at infinity is given by the limit 
u — > 0, i.e. X — > 00: 

u=~, y=-, z = - (5.5) 
In these coordinates, eqs. (5.1) take the form 

ii = -[(1 + gh)ez 2 + (1 + g 2 )p 2 ]u 7 

y = -[(l + gh)ez 2 + (1 + g 2 )p 2 ]y + (3 - h 2 )z 2 - (1 + gh)p 2 , (5.6) 
z = -[(1 + gh)ez 2 + (1 + g 2 )p 2 - y], 

where we have defined p = P/X and a dot denotes ud/d^. 

The discussion of the properties of the critical points at infinity is complicated, because 
it strongly depends on the range of variation of the parameters. The cases a = and 
g = —1/h are degenerate. In particular, the latter is completely integrable and is solved 
in the appendix. 

It is convenient to discuss separately the case of positive and negative A. 
A. A > 

For positive A, the critical points with u = can be classified in four categories: 
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1) A critical point, which we denote L, is placed at y = — ^j^p , z 2 = 0, with 
P 2 = TTP"- ^ ne eigenvalues of the system obtained linearizing (5.6) around this point are 

The point L is the endpoint of the trajectories lying in the Z = plane. The analysis of 
stability shows that it acts as an attractor on the trajectories coming from finite distance 
and, if 71 > 0, also on a two-dimensional bunch of trajectories lying on the surface at 
infinity (otherwise as a saddle point). 

2) If h 2 < 3, a critical point M lies at y = 77^, z 2 = (1 3 ~ g h fe 2 )2 , with p 2 = 0. The 
eigenvalues of the linearized system are — (double), — i+g h ■ 

This is the endpoint of the trajectories lying on the hyperboloid P = 0. The analysis 
of stability shows that, if 1 + gh > 0, it attracts the trajectories coming from finite distance 
and a one- or a two-dimensional bunch of those lying on the surface at infinity, depending 
on the sign of 72. If 1 + gh < 0, it repels all the trajectories. 



3) If a, 71 and 72 all have the same sign, a critical point N lies at y = 1, 



71 

a 



with p 2 = 22-. The eigenvalues of the linearized system are —1, — \{1 ± a/1 + 87172 /a). 
This critical point exists also if 72 has opposite sign, but is not acceptable physically, since 
corresponds to Q 2 < 0. 

The point N is the endpoint of the hyperboloid (3.2) in the X = Y plane. The 
eigenvalues containing the square root are real and have opposite sign. Hence N acts as 
an attractor on the trajectories coming from finite distance and as a saddle point on the 
trajectories at infinity. 

4) If a > 0, two critical points <5i,2 he at y = — ^-f^p^, z 2 = 0, with p 2 = 0. The 
eigenvalues of the linearized system are 0, 2, y. If h 2 > 3, the critical values of y both have 
the same sign (negative if 1 + gh > 0) , otherwise have opposite sign. 

The points Q are the endpoints of the trajectories with P = Z = 0. They act as centers 
on the trajectories coming from finite distance, while their nature for the trajectories at 
infinity depends on the sign of 3 — h 2 and 1 + gh. 



B. A < 

For negative A, the critical points at u = can again be divided in four categories: 

1) The critical point L is still present, and has the same properties as for positive A. 

2) If h 2 > 3, a critical point M lies at y = 77^77, z 2 = (1 ^~ ft 3 )2 , with p 2 = 0. The 
eigenvalues of the linearized system are the same as for the point M, but now the sign of 
the eigenvalues is different: if 1 + gh < 0, the point attracts the trajectories coming from 
finite distance and acts either as an attractor or as a saddle point on those lying on the 
surface at infinity, depending on the sign of 72. If 1 + gh > 0, it repels all the trajectories. 

3) If a > 0, 71 > and 72 < 0, a critical point N lies at y = 1, z 2 = — with 
p 2 = — . The eigenvalues of the linearized system are the same as for the point N, and all 
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take negative values. In analogy with N, this critical point exists also if 72 > 0, but the 
corresponding solutions are not physical. 

4) If a > 0, the two critical points Qi^ are still present, with the same properties as 
for A > 0. 

Asymptotic properties of the solutions 

As mentioned above, spatial infinity corresponds to the limit £ — > £o 5 where £0 is a 
finite constant. It is easy to see that for £ — > £0, the functions x an d V behave as 

e x He _ eo |iA>o e v^^_^vo/v 

where vq = (1 + gh)zQ + (1 + g 2 )pl, and the subscript indicates the value taken by the 
function at the critical point. From this one can deduce the behavior of the metric functions 
ending at the critical points. It is useful to write them in terms of the Schwarzschild 
coordinate r defined in (3.8) as J e 2iy+2p d^. Straightforward algebraic manipulations lead 
to the asymptotic behaviors listed in the table, where (3 is defined as in (3.23): 





e 2p 


e 2v 


e 2 ^ 


L 


r 2/(l+s") 


r -2/(l+g') 


r -2g/(l+g') 


M 


r 2/(l+h 2 ) 


r 2/(l+h 2 ) 


r 2h/{l+h 2 ) 


N 


r 2(g+h) 2 /p 


r 2(1+g 2 -h 2 )/p 


r 4(g+h)/p 



These coincide with the asymptotic behavior of the exact solutions found in sect. 3: more 
precisely, point L corresponds to case B, A = 0, point M to case A, Q = 0, and point N to 
case C. Of course, they are also identical to those obtained in ref. [1] in the planar limit. 

From the results listed above, we can deduce the global properties of the solutions, 
which depend on the values of the parameters A, g and h. The following picture of the 
phase space emerges: solutions with regular horizons start at points of the hyperbola (or 
ellipse if a < 0) of critical points at finite distance such that Xq = Yq = b/(g + h) 2 , b > 0. 

If a > and A > 0, black brane solutions exist with asymptotic behavior either L, M 
(if h 2 <3),or N (if 71, 72 >0). 

If a > and A < 0, black brane solutions exist with asymptotic behavior either L, M 
(if h 2 > 3), or N (if 71 > 0, 72 < 0). 

If a < regular black brane solutions may be still present. However, in this case the 
trajectories do not necessarily intersect the ellipse of critical points at finite distance, so it 
is not possible in general to predict the existence of solutions with regular horizons, that 
would have to start from points with Xq = Yq = b/(g + h) 2 . In particular: 

If a < and A > 0, black brane solutions may exist with asymptotic behavior either 
L or N (if 71,72 < 0). 

If a < and A < 0, black brane solutions may exist with asymptotic behavior either 
L or M (if h 2 > 3). 

The points L attract the trajectories lying in the plane Z = (i.e. those with vanishing 
potential) and a one-dimensional bunch of trajectories with Z ^ 0; the points M attract 
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the trajectories lying on the hyperboloid P = (i.e. those with vanishing charge) and, 
in some range of values of the parameters, also a one-dimensional bunch of trajectories 
with P 7^ 0; finally, for A > 0, the points N may attract a one-dimensional bunch of 
trajectories, depending on the values of the parameters g and h. Notice that the points Q 
are simply the endpoints of the curve of critical points at finite distance and do not attract 
any trajectory. 

Each solution depends on the two free parameters b (related to the mass) and Q 
(electric charge). In general, a third parameter, in addition to b and Q, may be necessary 
to parametrize the solutions. However, as emerges from the discussion of the exact special 
solution in the appendix, consistency of the thermo dynamical interpretation may require 
that this parameter be related to the electric charge. Also interesting is the possibility 
that solutions with different asymptotics exist for given values of the parameters A, g and 



6. Final remarks 

Although in the general case it is not possible to obtain analytic planar solutions of the 
EML model, we have discussed their global properties and classified all the possible regular 
black brane solutions in four dimensions. Several cases may arise, depending on the range 
of the parameters A, g and h that define the model, but we have been able to show that 
only three possibilities exist for the asymptotic behavior of the solutions: two of them are 
common with the limit of vanishing charge or potential, while one is characteristic of the 
general case. For given values of the parameters, solutions presenting different asymptotic 
behavior can exist. No asymptotically flat, anti-de Sitter or Lifshits solutions arise, except 
in the trivial case h = 0. 

The general analytic solution can be found for vanishing charge or potential, and for 
special values of the parameters g and h. In addition, also some special solutions can 
be found for generic values of the parameters. Some of these solutions had already been 
obtained in the literature [2], but not in full generality and in awkward coordinates. 

Acknowledgements 

I wish to thank Mariano Cadoni for some useful discussions. 

APPENDIX 

When g = — h -1 , the system can be completely integrated. This has been partially 
done in [2], but with a choice of coordinates that obscures the structure of the solutions. 
Our discussion illustrates the possibility that for a given value of h solutions with different 
asymptotic behavior exist, as it has been deduced in the general case from the study of 
the phase space. 

For g = — the system (2.5)-(2.6) diagonalizes, 



h. 



x " = (l + h- 2 )e 2 *, 



(A.l) 



and is solved by 



Ae 2 " = 




4b 2 K 2 e 2b Z 



(A.2) 



(3-/i 2 )(l-e 2 <) 2 ' 



(l + /i- 2 )(l-K 2 e 2 ^) 2 ' 
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and 

with a, 6, c and K integration constants. Substituting in the constraint (2.12), and requir- 
ing the existence of a regular horizon, one gets 

(3-/i 2 )(l+ft 2 ) 
" = » = (1 - fe y c - f- 4 '") 

One must now distinguish the generic solutions with K ^ 1 from the special solutions with 
K=l. 

Defining r = J e 2u+2p d^ « (1 - e 2a£)-(i+/i 2 )/(i-/i 2 ) ? the K = 1 solution can be put in 
the Schwarzschild form (3.6) with 

_ (l + h 2 ) 2 X 2{l+4h 2_ hi)/{l + h 2 )2 / _ ^\ 

(3-/i 2 ) I r<V' 



2 _ / (1 + fr 2 )Q 2 2 (l-/» 2 ) 3 /(l+/» 2 ) 2 p 20 _ r -4^(l-^) 2 /(l+ft 2 ) : 

1 l/ /i 2 (3-/i 2 )A ' 



(A5) 



where 5 = (3 — /i 2 )/(l + h 2 ) and is a free parameter. These solutions exist for h 2 < 3 if 
A is positive and for h 2 > 3 if A is negative. For positive \i, they represent a black brane 

with domain wall asymptotics, parametrized by mass density m = j^yj ^^g p-AQ 2 and 
charge density Q. 

If K 7^ 1, one obtains in the same way 

-^^(-^(-^"~. 

/ (l + /l 2 )Q 2 r 2/(l+^) , .^W) 



e 2 * = r 2h ^ 1+h ^ (l + ^) 



Z, n 2h/(l+h 2 ) 



These solutions present a further free parameter v ^ 0, that can be associated to the scalar 
charge. For \i, v > 0, they describe a 3-parameter family of black branes with domain wall 
asymptotic behavior different from that of the K = 1 solutions, which are recovered in the 
singular limit v — > oo. The curvature is finite at the point r = /U 1//<5 , that can be identified 
with a horizon, while it diverges at r = and r = — is, which are singularities of the metric. 
The interpretation of this solution is however not easy, in particular for what concerns 



its thermodynamics, because of the factor 1/ \Jv{y + fj,) in R 2 , that implies that there is no 
ground state corresponding to v = 0. A similar situation occurs for the neutral black branes 
discussed in [14]. In our case, the problem can be solved imposing that v{y + n) = uQ 2 , 
with a an arbitrary normalization, thus reducing the number of free parameters. While 
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for spherically symmetric solutions this condition is dictated by the normalization of the 
volume element, in the planar case it must be imposed by hand. Notice that this choice 
implies that \i and v must be positive, in order to avoid naked singularities. 

We remark that, while the solutions (A. 6) have been obtained in [2] in an equivalent 
but rather involved form, the solutions (A. 5) are new. 

The special cases h 2 = 1, 3 should be studied separately, since for these values of 
the parameter h some degeneracies appear in our calculations, but the generalization is 
straightforward, and we do not consider it in detail. 
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